Draft version October 30, 2012 

Preprint typeset using I^TgX style emulateapj v. 08/13/06 

Online-only material: color figures 



<N 

o 

<N 
-1— > 
o 

o 

(N 
O 

^: 



DIFFERENTIAL EMISSION MEASURE ANALYSIS OF MULTIPLE STRUCTURAL COMPONENTS OF 
CORONAL MASS EJECTIONS IN THE INNER CORONA 

X. CHENG^'2'^ J. Zhang^, S. H. SaarS h M. D. DlNG^'^ 
^ School of Astronomy and Space Science, Nanjing University, Nanjing 210093, China 
^ School of Physics, Astronomy and Computational Sciences, George Mason University, Fairfax, VA 22030, USA 
^ Key Laboratory for Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, Nanjing 210093, China and 
^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

Draft version October 30, 2012 

ABSTRACT 

In this paper, we study the temperature and density properties of multiple structural components 
of coronal mass ejections (CMEs) using differential emission measure (DEM) analysis. The DEM 
analysis is based on the six-passband EUV observations of solar corona from the Atmospheric Imaging 
Assembly onboard the Solar Dynamic Observatory. The structural components studied include the 
hot channel in the core region (presumably the magnetic flux rope of the CME), the bright loop-like 
leading front (LF), and coronal dimming in the wake of the CME. We find that the presumed flux 
rope has the highest average temperature (>8 MK) and density {r^l.{) xlO^ cm~^), resulting in an 
enhanced emission measure (EM) over a broad temperature range (3 < T(MK) < 20). On the other 
hand, the CME LF has a relatively cool temperature (~2 MK) and a narrow temperature distribution 
similar to the pre-eruption coronal temperature (1 < T(MK) < 3). The density in the LF, however, 
is increased by 2% to 32% compared with that of the pre-eruption corona, depending on the event 
and location. In coronal dimmings, the temperature is more broadly distributed (1 < T(MK) < 
4), but the density decreases by ^35% to ~40%. These observational results show that: (1) CME 
core regions are significantly heated, presumably through magnetic reconnection, (2) CME LFs are a 
consequence of compression of ambient plasma caused by the expansion of the CME core region, and 
(3) the dimmings are largely caused by the plasma rarefaction associated with the eruption. 
Subject headings: Sun: corona — Sun: coronal mass ejections (CMEs) 
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1. INTRODUCTION 

Coronal mass ejections (CMEs) are perhaps the most 
spectacular form of solar activity, which expel large quan- 
tities of plasma (order of ~10^^ to ~10^^ kg) at a speed 
of h undreds of km s~^ with the faste st ones over 3000 km 
s-i (|Yashiro et al.ll2004l : IChenll201lf ). Detailed kinematic 
analyses find that the acceleration of a CM E mainly oc- 
curs in the lower corona (e.g., < 3.0 Rq] IZhang et al.l 
l200lf ). Subsequently, it propagates into interplanetary 
space, probably tak i ng the form of a mag netic cloud 
(jBurlaga et al.l 119821 : [Klein fc Burla"gal Il982( ). A 



mag- 
netic cloud is able to produce severe geomagnetic dis- 
turbances if it intera cts with the Earth's magnetosphere 
(jCosling et al.lll993[ ). 

White-light coronagraph observations in the past 
decades have revealed that many CMEs display a char- 
acteristic three-part structure: a bright loop-like lead- 
ing front (LF), a dark cavi ty underneath, and an em - 
bedded bright compact core (jllling fc Hundhausenl ll983h. 
When a pre-CME structure lifts off from the associated 
source region, it can cause the expansion and successive 
stretching of the overlying magnetic field lines to form a 
CME. At the same time, the surrounding plasma accu- 
mulates at the CME fr ont, thus e nhancing the plasma 
density in the CME LF (|Cheng et al. 2011). In the mid- 
dle corona (e.g., 3-10 R©), densities in the LFs are usu- 
ally in the order of 10^ to 10^ cm~^, which represent 
~ 10-100 times enhancement over the background corona 



at th ose heights (jCiaravella et al.l I2003L 120051 : ISchwennI 
l2006f ). Temperatures of the LFs at 1. 5 R(^ have also 
been inferred, rangin g from 6.0 x 10^ K (ICiaravella et al.l 
I1997D to 2.0 X 10^ K (jBemporad et al.ll2007[ ). The bright 
cores of CMEs are usually believed to originate from the 
filament material (the cool but dense plasma suspended 
in the tenuous corona; Gopalswamv et al. 2006). Using 
ultraviolet spectral data, lAkmal et al.l (2001) estimated 
densities in the bright core ranging from 1.4 X 10^ to 7.0 
X 10^ cm~^ at 1.3 R©. Density quickly decreases with 
increasing height and vary from 1.3 x 10^ to 4.0 x 10^ 
cm-3 at 3.0 R© (Ra ymond fc Ciaravellall2004D . 

A flux rope structure, involving a set of twisted mag- 
netic field lines around a central axis, is often used to in- 
terpret the three-part structure of a CME; for instance, 
the dark cavity and the bright core of the CME corre- 
spond to the whole flux rope a nd the magnetic dips of the 
flux rope, respectively (e.g.. iLow fc HundhausenI 119951 : 
IChenI 119961 : ICibson et al.l 120061 : iRilev et al.ll2008h " Such 
helical flux rope conflguration has been reconstructed 
using nonlinear force-free fleld models b ased on pho- 
tosph e ric vector niagneto g ram data ( e.g. , ICanou etal 



20101 ). 



tospn eric vector magnetogram data ( e.g.. lUanou et al 
2009 1 : ICheng et al.l I2OIOI : iGuo et al.l I2OIOI : iJing eTal 
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Evidence for the existence of the flux ropes has been 
found in in-situ solar wind data, which often show a 
large an gle rotation of the magnetic fleld in rnagneti c 
clouds feurl aga et all 119821 : iKlein Burlagal Il982f ). 
Direct evidence of the flux rope comes in the existence 
of a conspicuous channel structure in the inner coron a 
before and during a solar eruption (|Zhang et al.ll2012l ). 
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This channel initially appears as a twisted and writhed 
sigmoidal structure in high temperature passbands, for 
example in 131 A at --10 MK and 94 A at ^ 6 MK seen by 
the A tmospheric Imaging Assembly (AIA; iLemen et alJ 
l2Qllf ) on board the Solar Dynamic Observatory (SDO). 
The channel eyolyes toward a semi-circular shape in 
the slow rise phase and then erupts upward rapidly 
in the impulsive acceleration phase, producing the 
front- cavity- core c omponents of the resulting CME 
(|Zhang et all l2Q12f ). The role that the hot channel 
plays in the eruption process appears similar to that 
of flu x ropes in the modeling and simula t ions of CMEs 
(e.g [Che^^' ^1996"; "Chen & Shibatal (2OOO!; "Lin & ForbesI 
2QQQ|; iTo rok & Kliem 2005; Kliem & Torok 2006; 
Aulanier et aL 2010; JFan Gibson! 120071 : iFanI I2OIOI : 
Qlmedo Zha ng 2010^ 



In addition to the three components, coronal dim- 
mings are another interesting phenomenon with close 
connect ions to CMEs. Dimmings can be obs erved in soft 
X-ray s (ISterhng fc Hudson 19971), in EUV (jZarro et all 
I1999I : [Tho mpson et ,aL 1998), and even in Ha pass- 
bands (pTang et alJ [2003[ ) . The commonly accepted 
physical explanation for coronal dimmings is that they 
represent a density drop in the inner corona result- 
ing from the plasma esc ape or depletion in the wake 
of a CME (Thompson e t all [1998: Harr ison fc LvonsI 
l2000f ). although plasma heating could also play a role, 
as hot plasma becomes less visible to the instru- 
ments sensitive primarily to lower temper atures (e.g., 
iRobbrecht fc Wand [20To[ I Cheng et al.ll201lf ). Assuming 
that EUV emission lines are optically thin and tempera- 
tures of dimmings do not change signiflcantly, iJin et all 
(|2009f ) estimated a density depletion of ^bO% at the early 
stage of dimmings. Later on, the intensity of dimming 
region gradually recovers, with several possible causes 
identified, including heating of c onfined plasma in coro- 
nal loops (|McIntosh et al.l [20071 ). interchange reconnec- 
tions between open mag netic field and small coronal 
loops (lAttrill et al.l 12008 ). or outfiows from the transi- 
tion region (jjin et al.ll2009f ). 

Previous studies have revealed the properties of CME 
structures and associated dimmings to a certain extent. 
Nevertheless, detailed information on the density and 
temperature properties of these structures is still lack- 
ing. Recently, differential emission measure (DEM) anal- 
ysis has been applied to diagnose the physical prop- 
erties of a CME. Using Hinode/ElS spectroscopic ob- 
servations, iLandi et al.r(j2Q10i) reconstructed the DEM 
distribution of the CME core and found the plasma 
in the CME core was heated slightly during the CME 
eruption. From the DEM n iaps derived by an a lyzings; 
four SOHO/EIT bandpasses. IZhukov fc Aucherd (j2004t 
found that the temperature of the EUV dimmings mainly 
varied between logT ~5.0 and logT ~6.5. Before and 
after the dimmings, the average DEM level decreased 
without a changs:e in the ov erall temperature distribution 
(also see.'Ti an et"al]l2012[ ). Moreover, the DEM method 
has also been used to determine the three-dimensional 
density and temper ature structures in the quiet Sun 
(|Vasquez et al.|[2^010f ) and to investigat e the temperature 
evolution in the post-fiare loop systems (jReeves fc Weberl 
[2OO9I ). 

In this paper, we apply the DEM method to the lat- 
est SDO /AIA data, which provides an opportunity for 



making a significant improvement in understanding CME 
structures. The detailed thermal properties of multiple 
CME components, including the fiux ropes (seen as hot 
channels in AIA observations), the bright LFs, and the 
dimmings are analyzed. Instrument and data reduction 
are presented in Section 2. In Section 3, we show the 
results, followed by discussion and conclusions in Section 
4. 

2. INSTRUMENT AND DATA REDUCTION 

2.1. Instrument 

The AIA on board SDO images the solar atmosphere 
through ten passbands almost simultaneously, with a 
temporal cadence of 12 s, a spatial resolution of 1.2'', 
and a field of view (FOV) of 1.3Rq. Six of the filters 
cover EUV lines formed at coronal temperatures at 131 
A (Fe VIH, Fe XX, Fe XXIH), 94 A (Fe XVIH), 335 A 
(Fe XVI), 211 A (Fe XIV), 193 A (FeXH, FeXXIV), and 
171 A (Fe IX), respectively. The temperature response 
functions of these passbands, as shown in Figure [H indi- 
cate an effective temperature coverage from 0.6 to 20 MK 
(O'Dwyer et al. 2010; Lemen et al. 2011,). During a so- 
lar eruption, the 131 A and 94 A passbands are sensitive 
to the hot plasma from eruption core regions, while the 
other passbands are bet ter at viewing the c o oler LFs and 
dimm ing regions (e.g., ICheng et al.l I201II IZhang et al.l 
l2012f ). The multi-passband broad-temperature capabil- 
ity of AIA makes it ideal for constructing DEM models 
of the distinct CME structures. 

2.2. Method 

The observed flux Fi for each passband can be deter- 
mined by: 



Ri{T) X DEM(T)dT. 



(1) 



where the Ri(T) is the temperature response func- 
tion of passband and DEM{T) denotes the plasma 
DEM in the corona. In this work, we use the 
"xrt_dem_iterative2.pro" routine in SSW package to com- 
pute the DEM. This code was origi nally designed for 
Hinode/XRT (X-r ay Telescope) data (jColub et al.ll2004l : 
IWeber et al.ll2004 ). and h ere is modified slig!:htly to wo rk 
with AIA data (see also ISchmelz et all l2010L l2011al lbl: 
IWinebarger et al.ll201lf ). For more details and tests of 
this method, see the Appendix. 

2.3. Data Analysis 

In this paper, we analyze three well-observed CME 
events, which occurred on 2010 November 03, 2011 
March 08, and 2011 March 07, respectively. We use 
the CME event on 2011 March 8 as shown in Figure 
[2] to illustrate our analysis process. First, we use the 
"aia_prep.pro" routine to process the AIA images in six 
EUV passbands to 1.5-level, which g uarantees a relative 
coalignment accuracy less than 0.6^^ (|Aschwanden et al.l 
[20TT[ ). Then, we outline three distinct regions (shown 
as boxes in Figure [2]) to compute the DEM. Regions a 
and b correspond to the fiare region and the quiet-Sun 
region, respectively, while a portion of the hot channel 
is selected in region c. In each region, the DN counts 
in each of the six passbands are normalized by the ex- 
posure time and spatially averaged over all pixels in the 
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region. We use these averaged count rates as the input of 
"xrt_dem_iterative2" routine to calculate the DEM curve. 

Figure [3{a) shows the DEM result for the flare region. 
The black solid curve indicates the best-fit DEM solution 
to the observed fluxes. In order to estimate DEM uncer- 
tainties, we compute 100 Monte Carlo (MC) realizations 
of the data. For each MC simulation, the observed flux Fi 
in each passband is perturbed by an amount of ^, which is 
randomly drawn from Gaussian distribution with a sigma 
equal to the uncertainty in the observed flux. The uncer- 
tainty is obtained by "aia_bp_estimate_error.pro" routine 
(Boerner 2012, private communication). The DEM code 
is then rerun for each of 100 MC realizations. The 100 
MC solutions thus represent 100 "equivalent" solutions 
of the original data within the noise on each channel. We 
use a blue rectangle, as show in Figure [3l^a), to represent 
the region surrounding the best-fit solution that contains 
50% of the MC solutions. The region consisting of two 
red rectangles and a blue rectangle covers 80% of the MC 
solutions. The region including all of colored rectangles 
contains 95% of the MC solutions. Thus, the upper and 
lower ends of these colored rectangles can be regarded 
as estimates of the uncertainties in the best-fit solution, 
indicating how well the DEM is determined at a given 
temperature bin. 

The flare region DEM shows a broad temperature dis- 
tribution from -0.8 MK to 20 MK (or 5.9 < logT < 7.3), 
indicating that plasmas with a wide range of tempera- 
tures are present in the compact flare region. Here, we 
introduce a useful parameter that characterizes the over- 
all temperature of the plasma, i.e., the DEM- weighted 
average temperature defined as: 



/ DEM{T) X TdT 
jDEM{T)dT 



(2) 



Using this definition, the average plasma temperature in 
the flare region is —9 MK, indicating that high temper- 
ature plasma dominates the emission of the flare region. 
Further, the errors of the DEM solutions are very small 
in the temperature range of 5.9 < logT < 7.3, indicating 
that the DEM is well constrained by the AIA data over 
most of the temperature range of the flare region. 

The DEM profiles for the selected quiet-Sun region and 
the fiux rope region are shown in Figure [3l^b) and (c), re- 
spectively. We can see that the emission in the quiet- Sun 
region is dominated by plasma with lower temperatures 
(5.9 < logT < 6.5), over which the DEM is weh con- 
strained; the quiet-Sun region has an average tempera- 
ture at —2 MK. The peak DEM of the quiet-Sun region is 
almost one order of magnitude lower than that of the fiare 
region. The fiux rope region seems to have two plasma 
components: the lower temperature component is similar 
in temperature distribution to that of the quiet region, 
but there is an additional, well-separated component of 
high temperature plasma. The average temperature of 
the fiux rope is —8 MK, indicating that it is somehow 
significantly heated during the eruption process. 

We also calculate the total emission measure (EM) us- 
ing: 

(3) 



EM = J DEM{T)dT. 



It is noteworthy that the total EM of the fiare region 
is —10^^ cm~^, almost two orders higher than that of 



the quiet-Sun region of 10^^ cm~^. Due to the contri- 
bution of its high temperature plasma, the total EM of 
the fiux rope is up to —10^^ cm~^, one order higher than 
that of the quiet- Sun region although lower than that of 
the fiare region. These total EM values, as well as the 
DEM- weighted average temperature are also indicated in 
Figure [3l Note that all above integrations are carried out 
over the same temperature range 5.9 < logT < 7.3. 

2.4. Uncertainties 

It is well known that DEM inversion is ill-posed and 
technically fraught with perils. On the one hand, errors 
in DEM inversion arise from the uncertainties in the re- 
sponse function Ri(T)^ including non- ionization equilib- 
rium effects, non-thermal populations of electrons, mod- 
ifications of dielectronic re combination ra tes ow ing to fi- 
nite d ensity plasmas (e.g.. ISummersl[T974 : Bad nell et al.l 
l2003f ). and even radiative transfer effects ()Judg?l2010f ). 
Moreover, the filling factor of the plasma is unknown, 
affect ing density d eterminations. Considering these ef- 
fects, iJudgd (|2010f ) estimated an uncertainty of 20% for 
Ri{T) although it is still a lower limit. 

On the other hand, errors in DEM inversion also orig- 
inate in the uncertainties in the background determina- 
tion, which is very importan t for DEM analysis (e.g., 
lAschwanden fc Boernerl [201lf ). In order to obtain the 
true DEM distribution inside the fiux rope, the emission 
from the background needs to be removed from the ob- 
served fiux. For fiux ropes, the background is determined 
from the nearby quiet- Sun regions (white boxes in Fig- 
ure |4]), which are close to and have the same heliocentric 
distance (to ensure a similar path along the line of sight 
in the corona) as the selected fiux rope sub-regions. We 
further inspect the effects of different backgrounds on 
our results and find that the DEM profiles do not change 
significantly, but the resulting parameters, e.g., average 
temperature, total EM, and density, vary by ~30%. For 
the CME LFs and dimming regions, however, we use 
the observed fiuxes to calculate the DEM directly with- 
out subtracting the fiux from a nearby region, since the 
emission in the LFs and dimmings comes from a large re- 
gion along the line of sight. In order to reveal the DEM 
changes, we compare the DEM in the same location be- 
fore and after the LFs and dimmings formation. 

Finally, we note that a cool line component was 
missing i n the old response function of AIA 94 A 
passband (iFoster fc Testal 120111 lAschwanden fc Boernerl 
[20111 : iSchmelz et al.l l2011bi r and there were also con- 
cerns abou^_tlie_j£cura^ of the 131 A response func- 
tion (jSchmelz et al.ll2011bD . We use the revised response 
functions (updated on 2012 January 30), in which the 
cool lines have been added to 94 and 131 A passbands, 
to avoid the issue of the missing cool line component. 

3. RESULTS 

3.1. DEM of CME Flux Ropes 

As previously noted, a CME fiux rope usually appears 
as an isolated channel structure in 131 A and 94 A images 
(jCheng et al.l 120111 : IZhang et al.ll2012D . Three such fiux 
ropes are shown in Figure |H They are viewed in different 
orientations: the first one is largely seen along its axis, 
while the other two are mostly viewed from the side. For 
each event, we select three different sub-regions along the 
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flux rope to calculate their DEMs (black boxes in Figure 
Hj); the corresponding background regions are shown by 
the white boxes in Figure HI The three sub-regions for the 
2010 November 03 CME flux rope are shown in the upper 
left panel of Figure HI one located close to the rim of the 
flux rope (region a), the other two close to the center of 
the flux rope (regions b and c). The resulting DEMs are 
shown in Figure [5l We find that at the rim of the flux 
rope, the plasma has an average temperature of T ~8.3 
MK, while at the center of the flux rope, the plasma 
temperature is T ~9.0 MK. The average temperature 
appears to decrease from the center to rim of the flux 
rope. But, the difference may not be taken seriously 
given the large uncertainty in these values. 

Based on the total EM of the flux rope, we can es- 
timate its density, assuming that the depth of the flux 
rope along the line of sight is approximately equal to its 
width, which can be measured directly (as shown by the 
black lines in Figure H]) . We calculate the density n in 
the flux rope using: 

n = ^EM/l. (4) 

where / is the depth (or width) of the flux rope. Note 
that, the filling factor in the density calculation is as- 
sumed to be 1. At the rim of the flux rope, with the 
estimated width of -40 Mm, the total EM of 2.1 xlO^^ 
cm~^ corresponds to a density of 7.2 xlO^ cm~^ (region 
a). Toward the flux rope center, the density increases to 
— 1.3 X 10^ cm~^, which is mainly due to the enhanced 
total EM (regions b and c). These results imply that the 
flux rope has not only a temperature structure but also a 
density structure. The maxima of both the temperature 
and the density occur near the flux rope center, with val- 
ues tending to decrease in the regions away from there. 
The average temperature, the density, total EM, as well 
as the width used in the calculation are also indicated in 
the left-upper corner of Figure [51 

Unlike 2010 November 03 event, the flux rope on 2011 
March 08 appears as a semi-circular tube with two foot- 
points fixed in the photosphere. The middle upper panel 
of Figure m shows the three small sub-regions we selected 
for DEM analysis. Due to a lower count rate at the top 
of the flux rope (region a), the DEM curve is poorly 
constrained in all but a few temperature bins. Never- 
theless, the derived parameters for this sub-region are 
T -10.9 MK, n -7.4 xlO^ cm"^ (Figure Ei; a)). In con- 
trast, the DEM curves in the two legs of the flux curves 
are well constrained, with the errors smaller (Figure [6l^b) 
and (c)). Most of the emission originates from high tem- 
perature plasma (6.6 < logT < 7.2; T -9.0 MK). The 
calculated densities are also higher (1.2-2.4 xlO^ cm~^). 

The flux rope on 2011 March 07 is similar to the 2011 
March 08 event in terms of orientation, but is associated 
with a filament. Most of cool filament material is located 
at the bottom of the flux rope structure, as shown in the 
right bottom panel of Figure HI From the DEM results 
of regions a and c (Figure [7j), we find that a significant 
amount of the emission is dominated by plasma with a 
cooler temperature (T <7.0 MK), likely due to the pres- 
ence of much cooler and denser filament material (0.6-1.7 
xlO^ cm~^). Region b is far from the filament region, 
thus is not contaminated by the cooler chromosphere ma- 
terial. It has a higher average temperature (T —9.0 MK) 
and a lower plasma density (5.3 xlO^ cm~^), and has 



more emission from the high temperature component of 
the plasma. Thus, we conclude that, except when mixed 
with a much cooler and denser filament component, a 
flux rope during its eruption is typically of a structure of 
hot plasma with an average temperature of —10.0 MK 
and a density of 1.0 xlO^ cm~^. 

The 2010 November 03 flux rope has exceptionally high 
counts throughout its eruption, thus making it an ideal 
case for studying the temperature evolution of the flux 
rope. We select the sub-regions including the maximum 
intensity as the flux rope centroid, which are indicated 
by the boxes in Figure [51 We derived the DEM distribu- 
tion every 12 seconds from 12:14 UT to 12:17 UT, and 
show the evolution of the average temperature in Figure 
[H We find that the flux rope centroid is further heated 
during the eruption; i.e., T increased from —8.0 MK to 
—10.0 MK, as the flux rope rose up and accelerated; the 
full kinemat i c evol ution of this event can be found in 
I Cheng et al.l (|201l[ ). Note that, just for the purposes of 
investigating the relative change of the temperature, the 
background fluxes were not subtracted from the observed 
fluxes before calculating the DEM curves (We implicitly 
assume the background is small and remains roughly con- 
stant.). 



3.2. DEM of CME Leading Fronts 

The eruption of the flux rope can push against the 
overlying magnetic field, whose expansion generates a 
compression front, observed here as the LF. A LF can 
be best seen in running or base difference images. Fig- 
ure [TOl shows the 171 A or 211 A base difference images 
of the three CME events. The base difference images 
are obtained through subtracting a pre-event image at a 
fixed time (the base) from the current images. The LF 
structure in the EUV images is very similar to those in 
coronagraph images. However, it is not clear whether the 
brightening front is caused by enhanced plasma density 
or an increase of temperature (or some combination) . We 
address this issue in this section using DEM analysis. 

The selected locations for the LFs of the three CMEs 
are shown in the boxes in Figure [TOl we calculate the 
DEM distribution for each region. In order to reveal 
the DEM changes after the LF formation, the DEM of 
the same location but for the pre-eruption state (pre- 
LF) is also calculated. The result of the 2010 November 
03 LF is shown in Figure [TOb). We find that the most 
significant (and well constrained) part of the DEM lies 
only in the lower temperature bins (6.0 < logT < 6.4). 
The DEM indicates that the LF plasma has no significant 
component with logT >6.5. 

One interesting result of these analyses is that the 
shapes of the DEMs do not vary greatly before and af- 
ter the LF appearance (Figure [TTTa) and (b)), indicating 
that the LF has a similar temperature distribution as the 
pre-eruption state. The average temperature for the LF 
and pre-LF are both —2.1 MK. This implies that the LF 
is not significantly affected by the strong heating process 
occurring near the core region during the eruption of the 
flux rope; almost all the released thermal energy, possi- 
bly via the process of magnetic reconnection, is confined 
to the region beneath the CME LF. 

We calculate the EM ratio R to study the variation of 
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the total EM: 

jDEMMT)dT 
J DEM^,,_MT)dT ^ ^ 

where DEMi^f and DEMp^e-LF represent the DEM of 
the LF and the pre-LF state (i.e., before the LP forms), 
respectively. The EM ratio value R is indicated in Figure 
[TTlfb). It is evident that the EM for the selected region 
increases up to ~13% when the LF appears. Since the T 
of LF and pre-LF are similar, this result strongly suggests 
that the brightening of the LF in EUV passbands is due 
to enhanced plasma density at the edge of the expanding 
CME rather than the change of temperature. In other 
words, the LF is a truly compression front, mainly the 
result of an enhancement in density. Assuming that the 
depth of the LF along the line of sight approximates its 
height from the solar surface (165 Mm), we find that 
the density increases from ^2.2 to 2.4 xlO^ cm~^. Note 
that, the estimated density is an upper limit, since the 
contribution from the background emission is included in 
the DEM calculation. 

The DEM distribution of the 2011 March 08 CME 
LF is similar to that of 2010 November 03 CME LF, 
with most emission coming from low temperatures: 6.0 
< logT < 6 A (Figure [TlTd)). Comparing with the pre- 
LF region, T of the LF does not change but the DEM 
temperature distribution broadens. Some emission also 
appears at low temperature (logT ~6.0) but this is less 
constrained (Figure [TlTc) and (d)). These changes result 
in the EM increasing by ^3% and the density increasing 
from ~1.0 to 1.1 xlO^ cm~^ when the LF appears. 

As for the 2011 March 07 CME LF, Figure [rTTf) shows 
its DEM result. The DEM is poorly constrained in the 
range of 6.0 < logT < 6.5. For the selected region, the 
average temperature slightly increases from T ~1.7 MK 
to T -1.9 MK following the LF formation. The EM is 
enhanced by —76% and the density is increased to —4.6 
xlO^ cm~^, compared to the pre-LF of —3.5 xlO^ cm~^ 
(Figure fTTT e) and (f)). It is worth mentioning that these 
values carry a considerable uncertainty, given the large 
errors in the DEM solutions. 

Note that, all the density changes that we derive above 
include an additional uncertainty from the simple estima- 
tion of the LF's depth. Moreover, the density change in 
the LF depends on sampled different regions, e.g., near 
the nose or at the fiank. Nevertheless, we find that the 
percentage of the density increase is always less than 50% 
in our study, thus providing an upper limit. 

3.3. DEM of CME Dimming Regions 

The dimming regions caused by 2010 November 03, 
2011 March 08, and 2011 March 07 CMEs are shown in 
Figure [121 To investigate the DEM of the dimming re- 
gions, we avoid areas that include hot plasma, e.g., the 
hot fiux rope areas in the images. The selected dim- 
ming sub-regions are indicated by the boxes in Figure 
[T2I We use the same method, comparing their DEM dis- 
tributions with the same region just before the dimming 
region appears (pre-D). The results for the three events 
are shown in Figure [T5l 

For the dimming on 2010 November 03 (Figure [TST b)). 
the well-determined portion of the DEM profile lies pri- 
marily in the range of 6.0 < logT < 6.6, and T -2.2 MK. 
The same region prior to the dimming has a very similar 



temperature, but the DEM peak following the dimming 
decreases by almost one order of magnitude. The to- 
tal EM decreases by 57% (Figure a) and (b)). Also 
assuming that the depth of the dimming along the line 
of sight is comparable to its height, we estimate that the 
density decreases from —3.5 to 2.3 xlO^ cm~^. Similarly, 
the DEM of 2011 March 08 dimming is centered around 
T —1.9 MK with temperatures spanning 6.1 < logT < 
6.4. The total EM decreases by —63% and the density 
decreases to —1.8 xlO^ cm~^ (Figure [T3lf c) and (d)). 
These results also hold for the 2011 March 07 dimming 
(Figure [T3lfe) and (f)), with the main DEM distribution 
spanning from logT —6.0 to logT —6.8, with T at —1.8 
MK, with the decrease of EM by 61% and of the density 
to —7.5 xlO'^ cm~^. Note that the DEM reconstruction 
here is rather poorly constrained; the errors vary by sev- 
eral orders of magnitude in most of temperature bins. 
Low counts are the cause of the large uncertainties in 
the inferred DEM. Nevertheless, we can conclude with 
high confidence that the dimming is mainly caused by 
the decrease of plasma density in the region. 

4. DISCUSSION AND CONCLUSIONS 

We summarize in Tabled] various properties of different 
CME structures, including average temperature, width 
of the DEM curve at the 10% of peak value (using a single 
Gauss fitting to the DEM curve), maximum DEM, total 
EM, and the density. The quantitative results from the 
DEM analysis further support our previous result based 
on a qualitative argument: a CME cons ists of a high ten i- 
perature fiux rope and a cooler LF (C heng et al.ll20l"ll ). 
By tracking the centroid of the fiux rope in the 2010 
November 03 event, we find that the DEM- weighted tem- 
perature of the fiux rope increases as it accelerates out- 
ward. Observations of the curve-in o f the fiux rope le^s , 
the shrinkage of the post-fiare loops (|Cheng et al 11201 If ), 
and the presence of high energ y hard X-ray sources sur- 
round ing the fiux rope (|Glesener et al.l I2OIII : iGuo et all 
l2012f ), taken together, argue that magnetic reconnect ion 
taking place in the current sheet underneath the fiux rope 
is responsible for heating the plasma inside the fiux rope 
to high temperature (making it visible at 131 A and 94 
A) . A similar flux rop e heating scenario was suggested 
by iLandi et al.l (|2010[ ) , who also attributed the heating 
to magnetic reconnection. 

We estimate the density of CME flux ropes during the 
eruption to be as high as '^1 x 10^ cm~ ^, similar to that of 
coronal loops in the lower corona (e.g. JUgarte-Urra et al.l 
120051 : lAschwanden et al.l 120081 ). This implies that the 
flux rope must originate in the core fleld structure of the 
associated active region, and that the high density of the 
entire structure is kept throughout the early evolution of 
the flux rope. 

The AIA data are also effective for reconstructing the 
DEM of CME LFs. Most of the LF plasma is conflned to 
a low-moderate temperature range (6.1 < logT < 6.5), 
with little emission at lower or higher temperatures and 
a DEM- weighted temperature in the range of T ~1.7 
MK to 2.1 MK. Since the plasma is prohibited from 
moving across magnetic fleld lines, the energy released 
by magnetic reconnection is difficult to transfer into the 
CME LF. Therefore, during the magnetic reconnection, 
the CME ffux rope is heated but the LF remains about 
the same temperature as the quiet coronal loops (e.g., T 
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-1.0-3.0 MK; iSchmelz et alllmol [20m]lbh . 

For the LF regions, the enhancement of the over ah 
EM (4%-76%) at typical corona temperature indicates 
that the brightening of the LFs is mainly due to an in- 
crease in plasma density. Assuming no depth change 
along the line of sight when the LF passes, we estimate 
a density increase of 2% to 32% for different LF regions 
at ^ 1. 2-1 ARcD. Using polar ized brightness Mauna Loa 
data. [Bemporad et al.l (|2007[ ) estimated that the density 
of CME LFs increases about 35% over the background 
coronal density at ^1.6Rq. Using AIA data and the 
DEM analysis, but assuming no te mperature change be- 
fore and after the LF appearance, iKozarev et al.l ()201ll ) 
found that the densities increase by —12% and —18% 
at — 1.3i?0 in two different bright fronts, respectively. 
Generally, our results are consistent with previous esti- 
mations for the density change of the CME LFs in the 
low corona. 

Dimming regions show changes in the opposite sense 
as the LFs. The well-constrained portion of dimming 
region DEM is in the range of 6.0 < logT < 6.5. DEM- 
weighted temperature varies from 1.7 MK to 2.2 MK, 
similar to that of the LFs. The DEM- weighted tem- 
perature doesn't change much before and after the dim- 
ming. However, the peak DEM decreases by as much 
as a factor of 10, which results in a total EM decrease 
of —60%. Similarly, assuming no depth change along 
the line of sight before and after the dimming, the de- 
creased EM corresponds to the depletion of the density of 
—40%. This shows that the dimming is mostly caused by 
the densi ty rarefaction or depletion in the lower corona 
falso see.lThompson et"aIlll998l:lHarrison fc Lvonsll2000l: 
Harrison et al.l 120031: IZhukov fc Aucherd 120041 : iJin et al.l 
20091 : iTian eF al. 201^. 

We summarize our main conclusions from the DEM 
analysis for the three distinct CME structures below. 

1. The plasma in flux ropes generates significant emis- 
sions over a broad temperature range of 6.5 < logT < 
7.3. For three flux ropes studied here, the densities vary 
from 0.5 xlO^ to 2.4 xlO^ cm"^ and the DEM- weighted 
average temperatures are all above 8 MK. In one case, 
the DEM- weighted temperature even increases to —10 



MK as the flux rope rises up, probably due to the con- 
tinuous magnetic reconnection. The presence of filament 
material within the magnetic dips of the flux rope can 
cause cooler apparent temperatures due to the mixture 
of hot and cold plasmas along the line-of-sight (e.g., 2011 
March 07 event as shown in Figure |4|). 

2. In three selected CME LF regions, the emission 
mostly comes from the cool plasma (6.1 < logT < 6.5), 
with the DEM- weighted average temperature unchanged 
by passage of the LF. Comparing to pre-LF regions, the 
density increases by 2% to 32%. We can thus conclude 
that the brightening of the LFs largely due to the plasma 
compression at the CME front rather than the increase 
of the temperature. 

3. Cool plasma (6.0 < logT < 6.5) also dominates the 
emission in the dimming regions, again with no change in 
the DEM- weighted average temperature. The decreased 
EM, and resulting dimmings, is thus due to density de- 
pletion in the low corona. The density reduction for three 
selected dimming regions is by —35% to —40%. 

In short, DEM analysis appears to be an important 
tool to diagnose the temperature and density proper- 
ties for various structural components during CME erup- 
tions. These basic parameters can provide valuable infor- 
mation to guide future CME modeling and simulations. 
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APPENDIX 

The ^^xrt_dem iterative2.pro" in the SSW package i s a DEM reconstruction routine, originally developed by Mark 
A. Weber fe.g.. iGolub eTatl |2QQ1 IWeber et al.ll2QQl . It uses a forward-fitting method, in which a DEM profile is 
guessed and then folded through the response of each passband to produce predicted fluxes. This process is iterated 
using Levenberg-Marquardt least-squares minimization, until the predicted fluxes are close to the observed ones. The 
DEM profile is interpolated using N-1 spline functions, representing the degrees of freedom for N different passband 
observations, which are directly manipulated by the well-known and much-tested "mpfit.pro" routine from Craig B. 
Markwardt^. 

To assess the ability of "xrt_dem_iterative2.pro" to reproduce DEMs, we have simulated AIA observations with sev- 
eral input DEM models, and then applied the forward- fitting routine to simulated AIA data to derive the best-fit DEM 
solutions. We then compared the reconstructed DEMs with input model DEMs, and explored uncertainties by comput- 
ing 100 different MC realizations (adding random noise within an uncertainty obtained by "aia_bp_estimate_error.pro" 
to the simulated observations) and fitting these as well. The fluctuations of 100 MC solutions measure the confidence in 
the reconstructed DEMs; typically, lower scatter in 100 MC solutions reflects smaller uncertainty in the reconstructed 
DEM solution. 

The model DEMs that we have tested include the single and d ouble isothermal DEMs, single and double Gaussian 
DEMs, and an active region DEM from CHIANTI (version 5.2.1: iDere et al.lll997l : ILandi et al.ll2006D . We also discuss 
some specific representative cases in more detail. 



^ http://cow.physics.wisc.edu/ craigm/idl/idl.html 
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In the single isothermal cases (Figure [14]), we find that in the temperature range of log Tq =5. 8-7. 2, the input 
temperature is successfully located, with only minor "spreading" into adjacent temperature bins. Since the input delta 
function DEM is fitted with splines, some EM "spillage" is to be expected, and the results are still quite good. The 
full width at half maximum (FWHM) of the recovered DEM in this To range is equal to 0.2 everywhere, but increases 
to 0.4 for log To out of this range. 

In the double isothermal cases (Figure [15]), we find that the fit EMs again always show some spread around the 
input temperatures logTi and logT2, and the FWHM tends to be larger (> 0.2). The two DEM peaks can not be 
resolved when the peak separation AlogT (logTi — logT2) < 0.4 (Figure [TSfb)). Nevertheless, the code still finds the 
DEM peak positions weh once AlogT > 0.4, independent of the EM1/EM2 ratio (Figure [ISj^a), (c), and (d)). The 
output parameters of fitting a double Gaussian profile to the double isothermal DEM for four specific cases are listed 
in Table [21 The recovered DEM curves can be always fitted successfully by the double Gaussian profiles with fixed 
log Ti and log T2 as same as modeled log Ti and log T2 except the case b in Figure [151 

In the single Gaussian DEM cases (Figure [16]). we find that for any log To in the range of 5.7-7.3 with an input 
width (Jt > 0.1, the best fitted DEMs are very close to the model Gaussian DEMs. The normalized error A = 
^{\DEMi - T)T;Mf^^^^^|)/ET)£;Mf'^^^^ < 0.12, the linear correlation coefficient CC > 0.99, and the chi-squared 
statistic = ^{DEMi- DEM^'''^''^)^ /a'^{DEMi) < 1.60 (see Table[3lfor metrics for the cases in Figure[l6]). We also 
note that the 100 MC solutions show little scatter around the best-fit solution, indicating the emission measures are fit 
with low uncertainty. Similar to the single Gaussian cases, with logT in the range of 5.7-7.3, any AlogT separation, 
and (Jt, the best fitted DEMs recover the double Gaussian DEMs profiles successfully (Figure [T7|) . The metrics of the 
goodness-of-fit for model DEMs are shown in Table [3l However, it is noted that the fit quality deteriorates when (JTi 
and/or = 0.1 (Figure [TTTc) and (d)). The normalized error A > 0.19, the linear correlation coefficient CC < 0.97, 
the increases to be 40.75 for the case d of Figure [TTl for example. 

Finally, the CHIANTI active region DEM "active_region_oso6.dem" is accurately reproduced across the entire tem- 
perature range with quite good accuracy (Figure [18]). The normalized error A = 0.12, the linear correlation coefficient 
CC = 0.99, and the chi-squared statistic = 4.37. Therefore, these tested cases suggest that "xrt_dem_iterative2.pro" 
is a reliable routine that can reconstruct the DEMs with AIA data well, with the caveats that it has more difficulty 
with isothermal plasmas, and with separating nearby EM peaks. This is likely due to its use of smooth s plines to define 
the DEM, which naturally have difficulty with extremely sharp features. It is worth mentioning that ISchmelz et al.l 
([2009a. b) compar ed ^xrt_dem_iteratiye2.prQ " with the Markov-Chain Monte Carlo (MCMC) based DEM reconstruc- 
tion algorithm bv iKashvap fc Drake! (| 19981 ). and found that the two methods are in good agreement. The MCMC 
method, since it searches space more thoroughly, generally finds slightly lower x^ minima, and "spikier" solutions, 
but at the cost of considerably more computing time. 
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Fig. 1. — AIA instrument temperature response curves for the six coronal passbands: 131 A, 94 A, 335 A, 211 A, 193 A, and 171 A. 
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Fig. 2. AIA 131 A, 94 A, 335 A, 211 A, 193 A, and 171 A images of the solar eruption on 2011 March 8. The sub-regions used for 
further analysis, indicated by boxes a, b, and c, are the flare region, quiet-Sun region, and flux rope region, respectively. 
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Fig. 3. — DEM curves for the flare region, quiet-Sun region, and flux rope region of 2011 March 8 CME, whose positions are shown in 
Figure [21 The black solid lines are the best-flt DEM distributions. The blue rectangle represents the region that contains 50% of the MC 
solutions. The two red rectangles, above and below the blue rectangle, and the blue rectangle compose the region that covers 80% of the 
MC solutions. All of colored rectangles form the region containing 95% of the MC solutions. Note that, to derive total EM, the DEM is 
integrated over the temperature range of logT=5.9 to logT=7.3. 
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Fig. 4. Top: AIA 131 A images of three CME flux ropes occurring on 2010 November 03, 2011 March 08, and 2011 March 07, 
respectively. The black boxes (a, b, and c) in each panel show the selected sub-regions used to reconstruct the DEM curves. The sub- 
regions shown by the white boxes indicate the locations where the background emissions are taken. Bottom: AIA 304 A images showing 
the CME flux rope-associated fllament on 2011 March 07. 
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Fig. 5. — DEM curves of the sub-region a, b, and c of 2010 November 03 flux rope (left upper panel of Figure 0). The black solid lines 
and the colored rectangles have the same meaning as in Figure [3] In order to derive the total EM and density n, the DEM is integrated 
over the temperature range of logT=6.0 to logT=7.2. 
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Fig. 7. — Same as in Figure [5] but for 2011 March 07 flux rope. 
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Fig. 8. — AIA 131 A images of the CME flux rope that occurred on 2010 November 03. The black boxes denote selected center regions 
of the flux rope during the eruption. 
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Fig. 9. — Temporal evolution of the DEM-weighted average temperature in the center of the flux rope for the 2010 November 03 CME 
as shown in Figure [U 
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Fig. 10. — AIA 211 A (left and right panels) or 171 A (middle panel) base-difference images of the 2010 November 03, 2011 March 08, 
and 2011 March 07 CME events. Their base images are taken at 12:00 UT, 03:30 UT, 19:30 UT, respectively. The black box in each panel 
shows the selected sub-region used to reconstruct the DEM curve of CME LPs. 
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Fig. 11. — DEM curves of three selected sub-regions in the LPs (shown by the black boxes in PigureO. The upper and bottom panels 
display the results for the pre-LP and LP, respectively. The black solid lines and the colored rectangles have the same meaning as in Pigure 
[31 To derive the total EM, EM ratio, and density n, the DEM is integrated over the temperature range of logT=6.0 to logT=6.5. 
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Fig. 12. — AIA 211 A (left and right panels) or 171 A (middle panels) base-difference images showing the dimmings associated with 
the three CMEs. Their base images are taken at 12:00 UT, 03:30 UT, 19:30 UT, respectively. The white box in each panel indicates the 
selected dimming region. 
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Fig. 13. — DEM curves of the three selected dimming sub-regions (shown by the white boxes in Figure^. The upper and bottom panels 
display the results for the pre-dimming and dimming regions, respectively. The black solid lines and the colored rectangles have the same 
meaning as in Figure [3] In order to derive the total EM, EM ratio, and density n, the DEM is integrated over the temperature range of 
logT=6.0 to logT=6.5. 
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Fig. 14. — Reconstruction of single isothermal DEMs with logT=6.0 (a), 6.5 (b), 5.7 (c), and 7.3 (d). The green solid lines show the 
model DEMs, the black solid lines display the best-fitted DEMs, and the colored rectangles have the same meaning as in Figure O 
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Fig. 15. — Reconstruction of double isothermal DEMs with with varying temperature separations and the DEM peaks. The green solid 
lines, the black solid lines, and the colored rectangles have the same meaning as in Figure fl4l 
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Fig. 16. — Reconstruction of single Gaussian DEMs with different logTo and ctt- The green sohd hues, the black solid lines, and the 
colored rectangles have the same meaning as in Figure [T4l 
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Fig. 17. — Reconstruction of double Gaussian DEMs with different temperature separations, DEM peaks, and (Tt- The green solid lines, 
the black solid lines, and the colored rectangles have the same meaning as in Figure [14] 
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Fig. 18. — Reconstruction of the active region DEM from CHIANTI file "active_region_oso6.dem" . The green sohd hnes, the black solid 
lines, and the colored rectangles have the same meaning as in Figure fl4] 
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TABLE 1 

Properties of CME flux ropes, LPs, and dimming regions. 



Event 


Region 


T 


(51ogT^ 


Maximum DEM 


Total EM 


Density 






(MK) 




(x 10^0 cm-s k-i) 


(x 10^'^ cm-^) 


(x 10^ cm-3) 


Plux Rope 




a 


8.3 


0.9 


3.6 


2.08 


0.72 


CMEl 


b 


9.1 


1.3 


8.8 


9.05 


1.22 




c 


9.0 


1.3 


10.2 


10.40 


1.32 




a 


10.9 


0.4 


2.6 


1.31 


0.74 


CME2 


b 


9.5 


0.9 


7.4 


4.74 


1.18 




c 


9.0 


1.3 


15.8 


13.70 


2.38 




a 


6.8 


0.9 


24.5 


12.00 


1.72 


CME3 


b 


9.0 


1.3 


1.2 


1.41 


0.53 




c 


4.5 


0.4 


3.9 


2.43 


0.59 


Leading Pront 


CMEl 




2.1 


0.9 


7.9 


0.94 


0.24 


CME2 




1.7 


0.9 


2.8 


0.26 


0.11 


CME3 




1.9 


0.9 


0.5 


0.06 


0.05 


Dimming Region 


CMEl 




2.2 


1.3 


3.9 


0.67 


0.23 


CME2 




1.9 


0.9 


3.8 


0.44 


0.18 


CME3 




1.8 


1.3 


0.7 


0.08 


0.07 



" Width of the DEM curve at 10% of the peak value based on a single Gaussian fit to the DEM curve. 



TABLE 2 

Parameters of fitting a double Gaussian profile to recovered double isothermal DEM. 



Case 


logTi 




logTa 






r\ 


a 


6.2 


0.2 


6.7 


0.2 


1.3 


0.5 


b 


6.3 


0.3 


6.5 


0.3 


1.4 


0.5 


c 


6.0 


0.2 


6.9 


0.2 


1.1 


0.9 


d 


6.0 


0.4 


6.9 


0.2 


1.1 


0.8 



" Ratio of the EM integrated with the recovered DEM to with the input DEM for the first isothermal peak. 
^ Ratio of the EM integrated with the recovered DEM to with the input DEM for the second isothermal peak. 



TABLE 3 

Metrics for goodness of fitted DEMs fPiGUREs fTMTs]) . 



Case 




CC^ 




Ri 


R2 






Single 


Gaussian 






a 


0.12 


0.99 


0.42 


1.12 




b 


0.12 


0.99 


1.14 


1.12 




c 


0.12 


0.99 


1.07 


1.12 




d 


0.12 


0.99 


1.60 


1.12 




Double Gaussian 


a 


0.15 


0.98 


5.72 


1.15 


1.05 


b 


0.14 


0.99 


8.20 


1.15 


1.07 


c 


0.20 


0.97 


12.01 


1.14 


0.99 


d 


0.19 


0.97 


40.75 


1.08 


0.99 


Real DEM 




0.12 


0.99 


4.37 


1.12 





« Normalized error A = j^DEM^odli ^• 

^ Linear Pearson correlation coefficient. 

^ Chi-squared statistic X = 2^ a'^(DEM ) 



